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This  paper  presents  the  results  of  a study  on  non- 
Gaussian  characteristics  and  statistical  prediction  of 
peaks,  troughs  and  peak-to-trough  excursions  of  coastal 
waves.  Through  the  statistical  analysis  of  more  than  500 
records  obtained  during  the  growing  stage  of  a storm,  it  is 
found  that  the  parameters  associated  with  skewness  and 
kurtosis  involved  in  the  non-Gaussian  probability 
distribution  are  significant  for  predicting  wave 
characteristics.  These  parameters  are  expressed  as  a 
function  of  water  depth  and  sea  severity.  The  limiting  sea 
severity  above  which  wind— generated  coastal  waves  are 
considered  to  be  non-Gaussian  is  obtained.  The  effect  of 
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water  depth  and  sea  severity  on  the  statistical 
characteristics  of  wave  period  are  clarified.  The 
probability  density  functions  applicable  to  peaks,  troughs, 
and  peak-to-trough  excursions  of  waves  in  finite  water 
depths  are  developed  based  on  the  non-Gaussian  probability 
distribution.  Numerical  computations  of  the  newly  developed 
probability  distribution  are  carried  out  at  various  water 
depths  for  various  sea  severities.  The  results  of 
computation  show  that  the  probability  density  functions  of 
peaks  and  troughs  both  agree  reasonably  with  the  histograms 
constructed  from  measured  data  obtained  during  the  ARSLCE 
Project. 
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CHAPTER  I 


INTRODUCTION 

The  profile  of  wind-generated  waves  observed  in  the 
area  of  finite  water  depth  is  significantly  different  from 
that  observed  in  deep  water,  and  the  difference  is 
pronounced  in  severe  seas.  Time  histories  of  the  former 
show  a definite  excess  of  high  crests  and  shallow  troughs  in 
contrast  to  waves  in  deep  water  as  demonstrated  in  Figure 
1.1. 

Analyses  of  wave  elevations  measured  in  finite  water 
depth  have  shown  that  the  histogram  of  wave  elevations  is 
not  symmetric  with  respect  to  the  mean  value,  and  its 
positive  skewness  increases  with  an  increase  in  sea  severity 
for  a given  water  depth.  This  implies  that,  from  the 
stochastic  point  of  view,  waves  in  finite  water  depth  can 
not  be  considered  a Gaussian  random  process  which  has  been 
proved  acceptable  for  probabilistic  prediction  of  waves  in 
deep  water. 

During  the  past  two  decades,  significant  progress  has 
been  made  in  the  probabilistic  prediction  of  various 
characteristics  of  random  waves  in  deep  water.  These  include 
prediction  of  wave  height,  period,  extreme  values,  breaking 
phenomenon  and  grouping  characteristics.  All  of  these 
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Example  of  wave  record  observed  in  finite  water  depth 


* 

Example  of  wave  record  observed  in  ocean 


Figure  1.1  Comparison  of  coastal  and  ocean  wave 
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prediction  methods  have  been  developed  based  on  the  Gaussian 
random  process  concept.  Needless  to  say,  this  concept  can 
not  be  applied  for  predicting  wave  characteristics  in  finite 
water  depth. 

Although  many  studies  have  been  made  on  wave 
characteristics  in  finite  water  depth,  only  a limited  number 
of  studies  has  been  carried  out  on  a method  for  predicting 
statistical  properties  of  wind-generated  waves  based  on  the 
non-Gaussian  random  process  concept.  In  particular,  the 
extent  to  which  coastal  waves  deviate  from  a Gaussian 
process  has  not  yet  been  clarified. 

The  purpose  of  this  s£udy  is  (1)  to  clarify  the  non- 
Gaussian  properties  of  coastal  waves  which  depend  on  water 
depth  as  well  as  sea  severity,  and  (2)  to  develop  a method 
for  predicting  wave  height  based  on  the  non-Gaussian  random 
process  concept. 

To  achieve  the  first  goal  of  this  study,  a 
statistical  analysis  is  carried  out  on  wave  time  histories 
obtained  by  the  Coastal  Engineering  Research  Center  (CERC) 
during  a storm  in  the  Atlantic  Ocean  Remote  Sensing  Land- 
Ocean  Experiment  (ARSLOE ) Project.  The  wave  time  histories 
were  obtained  at  ten  locations  in  water  depths  ranging  from 
1.35  m to  24.4  m along  the  CERC  Field  Research  Facility, 
Duck,  North  Carolina. 

For  the  second  phase  of  the  study,  probability  density 
functions  applicable  for  predicting  peaks  and  troughs  are 
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developed  separately,  since  the  magnitude  of  wave  peaks 
in  finite  water  depth  is  generally  much  greater  than  that  of 
troughs.  The  probability  density  function  of  peak-to-trough 
excursions  is  then  derived  assuming  that  peaks  and  troughs 
are  statistically  independent. 

This  paper  consists  of  three  chapters  followed  by 
conclusions  and  an  appendix.  Chapter  II  reviews  currently 
available  literature  on  probability  distribution  of  non- 
Guassian  random  processes  as  well  as  their  application  for 
predicting  statistical  characteristics  of  waves  in  finite 
water  depth. 

Chapter  III  discusses  the  results  of  analyzing 
measured  wave  data  based  on  the  non-Gaussian  probability 
density  function.  It  is  found  that  two  parameters  associated 
with  skewness  and  kurtosis  involved  in  the  non-Gaussian 
probability  density  function  are  significant  for  predicting 
wave  characteristics  in  finite  water  depth.  These  two 
parameters  are  expressed  as  a function  of  sea  severity 
and  water  depth.  The  limiting  sea  severity  below  which  the 
wind-generated  coastal  waves  are  considered  to  be  Gaussian 
is  obtained  for  a given  water  depth.  The  effect  of  water 
depth  and  sea  severity  on  the  statistical  characteristics  of 
wave  period  is  clarified. 

Chapter  IV  presents  the  development  of  a method  for 
predicting  peaks  and  troughs  as  well  as  peak-to-trough 
excursions  of  waves  in  finite  water  depth  based  on  the  non- 
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Gaussian  probabiltiy  density  function.  Numerical 
computations  are  carried  out  at  water  depths  ranging  from 
1.35  m to  8.77  m for  various  sea  severities.  The  results  are 
compared  with  measured  data  obtained  during  the  ARSLOE 
Project. 

Chapter  V summarizes  the  findings  and  conclusions 
obtained  from  the  results  of  the  present  study. 

In  Appendix,  the  wave  data  used  in  the  present  study 
are  described. 


CHAPTER  II 


REVIEW  OF  LITERATURE 

This  chapter  reviews  currently  available  literatures 
on  the  probability  distribution  applicable  to  non-Gaussian 
random  waves  as  a prelude  to  the  present  study. 

Ochi  (1985)  has  recently  presented  a state-of-the-art 
review  on  stochastic  analysis  and  probabilistic  prediction 
of  non-Gaussian  random  processes  in  ocean  engineering  in 
which  he  grouped  the  derivation  of  the  probability  density 
functions  applicable  to  non-Gaussian  random  process  into  two 
categories:  one  derived  from  probability  theory,  the  other 
derived  from  applying  Stokes  wave  theory.  In  particular,  the 
latter  approach  is  concerned  with  wave  surface  profiles 
which  are  expressed  as  an  accumulation  of  profiles  expressed 
in  the  form  of  a series  rather  than  accumulation  of  simple 
sinusoidal  form  on  the  wave  profile.  Hence,  this  approach 
actually  imposes  a preliminary  form  on  wave  surface 
profiles . 

There  are  three  probability  distributions  developed  by 
applying  the  concept  of  probability  theory.  These  include 
the  Gram-Charlier,  Edgeworth,  and  Longuet-Higgins 
probability  density  functions.  Gram-Charlier  probability 
density  function  ( Cramer , 1 946 ; Kendal  and  Stuart, 1961)  is 
derived  from  the  concept  of  orthogonal  polynomials.  Since 
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Hermite  polynomials  are  orthogonal  with  respect  to  the 
normal  distribution,  the  probability  density  function  is 
expanded  in  a series  of  Hermite  polynomials.  The  Edgeworth 
series  probability  density  function  (Edgeworth  1905)  is 
derived  following  the  same  concept  as  the  Gram-Charlier 
series  but  the  series  expansion  method  is  different  and  the 
expanded  terms  are  expressed  in  term  of  cumulants.  On  the 
other  hand,  the  Longuet-PIiggins  series  probability  density 
function  ( Longuet-Higgins  1963)  is  derived  by  applying  the 
cumulants  generating  function.  In  his  approach,  the 
integrations  involved  in  the  derivation  are  expressed  in 
terms  of  Hermite  polynomials. 

It  is  noted  that  these  three  probability  distribution 
functions  can  be  reduced  to  the  same  form  shown  below, 
although  their  method  of  derivation  are  entirely  different. 


f (x)  = 


/ 2 ir 


2a- 


1 + 


3 1 


V — ’ + 


4 i 


H ( — ) 

4 a 


5 1 


5 X 

- V— > + 1 


6 i 


2 
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) Ks(  — ) 


+ . 


(2.1) 
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where 

02  = variance  of  x 

_ 3 

A 3=  K3/(/k  ) = skewness 

_ 4 

^4=  K /( /k2 ) = kurtosis  - 3 

_ 5 

a5=  k5/(/k2) 

K . = cumulants 

3 

. X 

H^(— — ) = Hermite  polynomial  of  degree  n 


The  first  term  in  the  series  is  the  normal  distribution;  all 
other  terms  in  the  series  are  attributed  to  the  effect  of 
non-Gaussian  characteristics. 

Another  approach  for  deriving  the  probability 
distributions  applicable  to  non-Gaussian  random  process  is 
based  on  the  concept  that  the  wave  surface  profile  can  be 
represented  approximately  by  the  Stokes  wave  given  by 


1 1 3 

n = a2k  + a cosX  + a2k  cos2X  + a3k2cos3X 

2 2 8 

+ (2.2) 

where 


n = wave  surface  elevation,  a = a(  x,  t)  = amplitude 
X = kx  - wt  + e e = random  phase 

k = w2 / g = wave  number 


Tayfun  (1980)  developed  the  probability  density 

function  of  wave  surface  elevations,  n,  based  on  Stokes 

wave  2nd  order  expansion  in  deep  water  without  the  constant 
l 2 

term,  _ a k,  which  is  the  second  and  third  term  in 
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Equation  (2.2).  By  assuming  that  the  amplitude,  a,  follows 
Rayleigh  distribution  and  using  transformation  of  random 
variables,  the  cumulative  distribution  function  of  surface 
wave  elevation,  F(n  ),  is  derived.  A numerical  computation 
is  required  to  evaluate  the  probability  density 

function,  f (n  ) . 

Huang  et  al.(1983)  applied  Stokes  wave  of  3rd  order  in 
deep  water  which  is  the  first  four  terms  in  Equation  (2.2) 
to  represent  the  wave  surface  profile.  The  probability 

density  function  of  wave  surface  elevations,  n , is  derived 
by  assuming  the  steepness  of  wave  is  very  small  and  then 
using  a perturbation  scheme.  Huang  et  al.  (1983)  further 
derived  the  probability  density  function  of  n in  finite 
water  depth.  For  this,  the  Stokes  waves  of  3rd  order  in 

finite  water  depth  is  considered.  It  is  given  by  Bowden 

(1948)  as  follows: 
a2  k 

n = coth  Kd  + a cos  X 

2 

a2  k 3 

+ coth  kd  ( 1 + ) cos2X 

2 2 sinh2kd 

3 3 2 3 3 9 

+ a k ( 1 + + + ) cos3X 

8 sinh2  kd  sinh3 4 *kd  8 sinh6  kd 

(2.3) 

where 

d = water  depth 
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The  probability  density  functions  of  amplitudes  (peaks 
and  troughs)  and  wave  heights  have  been  derived  from  the 
same  two  approaches  as  considered  for  the  derivation  of  the 
probability  density  function  of  wave  surface  elevations: 
one  from  probaility  theory,  the  other  from  applying  Stokes 
wave  theory. 

Bitner  (1980)  derived  the  probability  density  function 
of  amplitude  and  wave  height  by  assuming  that  the  wave 
surface  elevation,  n , can  be  expressed  as 

n = a cos  ( wt  - e(t))  (2.4) 

where 

a = a(t)  = random  amplitude,  e(t)  = random  phase 

Then,  by  applying  Equation  (2.1)  for  the  probability  density 
function  of  n,  the  probability  density  function  of 

amplitude  is  derived.  Unfortunately,  the  probability  density 
function  of  amplitudes  developed  by  Bitner  is  not  applicable 
for  waves  in  finite  wate  depth,  since  the  distributions  of 
peaks  and  troughs  are  substantially  different,  in  general, 
as  demostrated  in  an  example  shown  in  Figure  1.1.  It  is 
highly  desirable  to  derive  the  probability  density  functions 
applicable  for  peaks  and  troughs  separately  for  non-Gaussian 
random  waves. 

Arhan  and  Plaisted  (1981)  developed  the  probability 
density  functions  of  peaks  and  troughs  based  on 


the 


} 


11 

application  of  the  Stokes  wave  of  2nd  order  in  finite  water 
depth.  They  expressed  the  peaks  and  troughs  of  Stokes  wave 
as  follows: 


£ = a + eB  a1 2 

where 

£ = +1  for  peak,  -1  for  trough, 

k 3 

B = coth  kd  ( 1 + ) 

2 2 sinh2kd 

d = water  depth 


(2.5) 


Here,  the  amplitude,  a,  follows  the  Rayleigh  distribution. 

Tayfun  (1984)  developed  the  probability  density 
functions  of  peaks  and  troughs  by  applying  the  Stokes  wave 
of  2nd  order  in  deep  water  with  a parameter,  a , controlling 
the  effect  of  the  nonlinear  component.  That  is. 


1 

n = a cosX  + aka2cos2X  (2.6) 

2 


where 

n = wave  surface  elevation 


0 £ a £ i 
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The  peaks  and  troughs  associated  with  the  wave  profile  shown 
in  Equation  (2.6)  can  be  expressed  as 


where 

e 


1 

£ = e a + a ka2 

2 


!+  1 for  peaks 

- 1 for  troughs 


(2.7) 


The  amplitude f at  follows  • the  probability  distribution 
developed  by  Tayfun  (1981)  which  considered  the  breaking- 
wave  effect. 

Tayfun  (1983)  also  developed  the  probability  density 
function  of  peak-to-trough  excursions  of  non-Gaussian  waves 
based  on  the  application  of  Stokes  wave  of  2nd  order  in  deep 
water  which  is  the  second  and  third  terms  in  Equation  (2.2). 

As  this  brief  review  of  the  literature  indicates, 
several  studies  have  been  carried  out  for  the  developement 
of  probability  density  functions  of  wave  surface  elevations, 
peaks,  troughs,  and  peak-to-trough  excursions  based  on  the 
Stokes  wave  theory  in  representing  wave  surface  profile. 
This  is  certainly  a way  to  derive  the  probability 
distribution  applicable  to  non-Gaussian  random  waves. 
However,  this  approach  imposes  a preliminary  form  as  a 
function  of  frequency  on  wave  profile  which  results  in  a 
limitation  on  the  non-Gaussian  characteristics  of  random 
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waves  in  finite  water  depths.  It  seems  to  be  desirable  to 
develop  the  probability  density  functions  of  amplitude 
(peaks  or  troughs)  and  wave  heights  by  applying  the  concept 
of  probability  theory,  which  does  not  impose  any  preliminary 
form  on  the  wave  surface  profile. 


CHAPTER  III 


NON-GAUSS IAN  CHARACTERISTICS  OF  COASTAL  WAVES 

As  stated  in  the  introduction,  the  profile  of  wind- 
generated waves  observed  in  coastal  water  shows,  in  general, 
a definite  excess  of  high  crests  and  shallow  troughs  which 
is  significantly  different  from  waves  observed  in  deep 
water.  This  difference  is  particularly  pronounced  in  severe 
seas.  This  implies  that  waves  in  coastal  water  where  the 
effect  of  water  depth  on  wave  characteristics  is  present 
cannot  be  considered  a Gaussian  random  process.  Hence, 
clarification  of  the  non-Gaussian  characteristics  of  coastal 
waves  is  desirable  to  provide  a basis  for  statistical 
prediction  of  wave  characteristics. 

The  purpose  here  is  to  clarify  the  non-Gaussian 
properties  of  coastal  waves  which  depend  on  water  depth  as 
well  as  sea  severity.  This  chapter  is  in  four  sections.  The 
first  section  summarizes  results  of  the  statistical 
analysis  based  on  the  non-Gaussian  probability  density 
function.  From  analysis  of  more  than  500  wave  records 
measured  at  various  water  depths  during  the  growing  stage  of 
the  storm,  it  is  found  that  two  parameters  associated  with 
skewness  and  kurtosis,  respectively,  involved  in  the 
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Longuet-Higgins ' non-Gaussian  probability  density  function 
are  significant  for  predicting  wave  characteristics  in 
finite  water  depth.  The  second  section  presents  the  analysis 
results  in  which  these  two  parameters  are  expressed  as  a 
function  of  water  depth  and  sea  severity.  The  third  section 
discusses  the  limiting  sea  severity  below  which  wind- 
generated coastal  waves  can  be  considered  as  a Gaussian 
random  process.  The  fourth  section  discusses  the  statistical 
characteristics  of  the  period  of  coastal  waves. 


3 . 1 Probability  Distribution  of  Wave  Profiles 

Several  probability  density  functions  have  been 
developed  for  representing  non-Gaussian  random  processes  as 
outlined  in  the  previous  chapter.  Since  an  important 
objective  of  the  present  study  is  to  clarify  the  extent  to 
which  coastal  waves  deviate  from  a Gaussian  random  process, 
it  may  be  well  not  to  impose  any  preliminary  form  on  wave 
profile  representation  such  as  Stokes  expansion  in 
analyzing  data  obtained  in  coastal  water.  Therefore,  the 
following  non-Gaussian  probability  density  function  derived 
by  Longuet-Higgins  is  considered  in  the  statistical  analysis 
of  wave  profiles. 
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a2=  variance  of  wave  deviation  from  the  mean 

a3=  k3  /(/ k2 )3 

A 4 = K4  /(/K2)4 

. *5=  % /(/-K2)5 

K.  = cumulants 
D 


x 

) = Hermite  polynomial  of  degree  n 


It  is  noted  A3  is  the  skewness  of  the  random 
variable  x,  and  a4  is  equal  to  the  kurtosis  minus  3.  The 
value  of  A3  and  A4may  be  positive  as  well  as  negative  in 
general;  however,  positive  values  are  observed  for  most  wave 
profiles  in  the  coastal  water. 

In  applying  the  above  probability  density  function  for 
the  analysis  of  coastal  waves,  the  following  remarks  are 
made : 

(1)  In  general,  the  accuracy  of  a function  which  is 
expressed  in  the  form  of  a series  of  uniform  convergence 
increases  with  increase  in  higher  order  terms.  However, 
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this  is  not  the  case  for  the  probability  density  function 
given  in  Equation  (3.1).  The  results  of  comparisons  between 
histograms  and  the  probability  function  have  shown  that 
higher  order  terms  do  not  necessarily  yield  better 
agreement.  Therefore,  it  is  desirable  to  examine  terms  of 
the  series  which  significantly  contribute  to  the 
distribution. 

(2)  The  probability  density  function  given  in  Equation 
(3.1)  at  times  becomes  negative  for  large  negative  x 
depending  on  the  value  of  a2  and  the  number  of  terms  taken 
into  consideration  from  the  series.  However,  this  will  not 
cause  any  trouble  in  practice,  since  the  x value  where  the 
function  becomes  negative  is  usually  outside  the  range 
where  the  histogram  exists  as  will  be  shown  later. 

In  order  to  clarify  possible  questions  arising  from 
the  above  remarks,  more  than  500  wave  records  measured  at 
various  water  depths  during  the  growing  stage  of  the  storm 
during  the  ARSLOE  project  are  analyzed.  Details  of  wave  data 
used  in  the  present  study  are  summarized  in  the  Appendix. 

Figure  3.1  shows,  as  an  example,  a comparison  between 
the  probability  density  function  given  in  Equation ( 3 . 1 ) with 
a histogram  of  the  wave  profile  obtained  at  an  average  water 
depth  of  1.35  m.  The  significant  wave  height  is  2.10  m 
which  is  considered  severe  for  this  water  depth.  Four 
variations  of  the  probability  density  are  compared  with  the 
histogram.  These  four  variations  represent,  respectively, 
the  first  term  of  the  series  (which  is  the 


normal 


PROBABILITY  DENSITY  IN  I/M 


1.2 


Figure  3.1  Comparison  between  observed  histogram,  Gaussian 
distribution  (heavy  line) , and  non-Gaussian 
distribution : 

GAGE  615,  average  water  depth  1.35m, 

water  depth  2.28  m at  13:15  on  Oct.  25th  due  to 

tide,  significant  wave  height  2.10  m 
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distribution),  the  first  two  terms  (including  the  parameter 

\3  ),  the  first  three  terms  (including  the  parameters  I3  and 
2 

A3  ),  and  the  first  four  terms  (including  the  parameters 
2 

A 3 , A 3 , and  A 4 ) . 

As  can  be  seen  in  this  figure, 

(1)  The  histogram  deviates  substantially  from  the 
normal  distribution. 

(2)  The  probability  density  function  becomes  negative 
for  large  negative  x.  However,  the  magnitude  of  the  negative 
probability  density  is  relatively  small,  on  the  order  of  2 
percent  or  less.  Furthermore,  the  negative  density  occurs 
outside  the  range  of  the  histogram.  Therefore,  it  will  not 
cause  any  serious  problem  if  we  assume  this  negative 
probability  density  to  be  zero  and  in  turn  normalize  the 
entire  probability  density  function  so  that  the  area  of  the 
density  function  becomes  unity. 

(3)  The  probability  density  function  which  includes 
the  parameter  A3  only  agrees  reasonably  well  with  the 
histogram  and  differs  substantially  from  the  normal 
distribution. 

(4)  The  agreement  with  the  histogram  becomes  poor  if 


2 


the 

term  with 

the  parameter 

A 3 

is  included 

in 

addition 

to 

the 

A 3 -term 

in  the  probability 

density  function.  Although 

not 

included 

in  Figure  3.1, 

the 

: same  trend 

is 

obtained 

if 

the 

term  with 

the  parameter 

a4 

is  included 

in 

addition 

to 

the  A 3 -term. 
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(5)  The  probability  density  function  consisting  of  the 

. 2 
first  three  terms  taking  the  parameters  X3,  X3,  and  X4 

into  consideration  agrees  well  with  the  histogram. 

(6)  The  addition  of  the  parameter  ^5  to  the 

probability  density  function  does  not  yield  any  appreciable 
change  in  the  shape  of  the  probability  density  function. 

The  same  trends  stated  above  was  also  observed  for 

many  other  wave  records  obtained  at  various  water  depths  in 

various  sea  severities.  Some  examples  are  shown  in  Figure 
3.2(a)  through  3.2(g)  for  which  the  water  depth  ranges  from 
1.35  m to  8.77  m.  In  each  figure,  five  variations  of 

probability  density  function  are  compared  with  the 
histogram.  These  five  variations  represent  different 
combinations  of  terms  in  Equation  (3.1).  They  are 

(1)  the  first  term  (GAUSS  in  the  figure), 

(2)  the  first  two  terms  ( NG-1 , triangle  marks  in  the 
figure ) , 

(3)  the  first  three  terms  (NG-2,  cross  marks  in  the 
figure) , 

(4)  the  first  four  terms  (NG-3,  star  marks  in  the 
figure) , 

(5)  the  first  five  terms  (NG-4,  square  marks  in  the 
figure ) . 

From  the  results  shown  in  the  figures,  we  may  safely 
conclude  that  the  parameter  X3  which  represents  the 
skewness  of  the  wave  profile  is  the  dominant  parameter 
affecting  the  non-Gaussian  characteristics  of  waves  in 


C X ) 
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Figure  3.2  Comparison  between  observed  histogram,  Gaussian 

distribution  (heavy  line),  and  non-Gaussian  distribution: 
(a)  Average  water  depth  1.35  m,  significant  wave 
height  1.60  m. 
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Figure  3.2  (continued) 

(b)  Average  water  depth  1.97  m,  significant  wave 
height  1.89  m • 
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Figure  3.2  (continued) 

(c)  Average  water  depth  3.70  m,  significant  wave 
height  2.44  m 


•CX) 
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Figure  3.2  (continued) 

(d)  Average  water  depth  5.55m,  significant  wave  height 
2.09m  ■' 


( X ) 
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Figure  3.2  (continued) 

(e)  Average  water  depth  5.97  m,  significant  wave  height 
2.18  m 
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Figure  3.2  (continued) 

(f)  Average  water  depth  6.97  m,  significant  wave 
height  3.17  m 


( X J 


27 


Figure  3.2  (continued) 

(g)  Average  water  depth  8.77  m,  significant  wave  height 
3.75  m 
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coastal  water  and  that  the  combination  of  the  parameter 
X3,  X 3 , and  X best  represents  the  non-Gaussian 
probability  density  function  given  in  Equation ( 3 . 1 ) . This 
conclusion  agrees  with  that  made  by  Huang  and  Long  (1980) 
in  their  experimental  study  of  surface  elevation  probability 
distribution.  On  the  other  hand,  the  combination  of  the 
parameters  A3  and  A3  or  the  combination  of  the  parameters 
A 3 and  a4  does  not  yield  satisfactory  agreement  between  the 
histogram  and  non-Gaussian  probability  density  function. 

Based  on  the  above-mentioned  conclusions,  only  the  two 
parameters  a3  and  x4  involved  in  the  non-Gaussian 
probability  density  function  are  considered  in  the  following 
statistical  analysis. 


3 . 2 Analysis  of  Parameters  of  the  Probability  Distribution 

It  was  shown  in  the  previous  section  that  X3  and  A4  in 
Equation  (3.1)  are  the  significant  parameters  which  govern 
the  non-Gaussian  probability  density  function  for 
representing  the  distribution  of  wave  surface  elevations. 
Since  the  non-Gaussian  characteristics  of  coastal  waves 
depend  on  water  depth  and  sea  severity,  it  may  be  well  to 
examine  the  properties  of  these  parameters  as  a function  of 
water  depth  and  sea  severity. 

It  was  first  thought  that  the  parameter  A3  (skewness) 
was  a function  of  significant  wave  height  and  wave  period. 
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However,  wave  records  show  that  X3  appears  to  depend  on 
significant  wave  height  only.  To  substantiate  this,  the  time 
variations  of  significant  wave  height  Hg,  parameter  X 3 wave 
period  (average  zero-crossing  wave  period)  and  water  depth 
obtained  at  average  water  depths  of  1.35  m and  8.77  m 
during  the  growing  stage  of  a storm  are  shown  in  Figure 
3.3(a)  and  3.3(b)  respectively. 

As  can  be  seen  in  Figure  3.3(a),  the  water  depth 
varies  with  time  by  a substantial  amount  due  to  the  tide. 
This  results  in  an  increase  and  decrease  in  significant  wave 
height  with  the  same  period  as  the  tide  (12  hrs.)  but  the 
magnitude  of  significant  wave  height  at  high  tide  increases 
consistantly  during  the  growing  stage  of  the  storm.  The  time 
history  of  the  magnitude  of  the  parameter  X3  (skewness) 
computed  at  hourly  intervals  (with  a few  points  missing  due 
to  bad  measurement)  demonstrates  almost  the  same  pattern  as 
that  of  the  significant  wave  height.  However,  the  magnitude 
of  the  zero-crossing  average  wave  period  slowly  increases 
with  time  during  the  growing  stage  of  the  storm. 

Another  example  observed  at  the  average  water  depth 
of  8.77  m is  shown  in  Figure  3.3(b).  The  example  shows  the 
water  depth  still  varies  with  time  by  a substantial  amount 
due  to  tide,  which  is  the  same  pattern  observed  in  the 
previous  example.  However,  the  variation  of  water  depth  with 
time  does  not  result  in  any  significant  increase  or  decrease 
in  the  magnitude  of  significant  wave  height  as  it  did  in  the 
previous  example.  The  magnitude  of  significant  wave  height 
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Figure  3.3  Variation  of  water  depth,  wave  period, 
significant  wave  height,  and  skewness 
with  time: 

(a)  Average  water  depth  1.35  m 
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Figure  3.3  (continued) 

(b)  Average  water  depth  8.77  m 
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increases  consistently  during  the  growing  stage  of  the 
storm.  A similar  trend  is  also  observed  in  the  time 
variations  of  parameter  ^ as  well  as  average  zero-crossing 
wave  period. 

From  the  results  of  analysis,  several  trends  are 
observed.  They  are 

(1)  The  time  variation  of  water  depth  due  to  tide 
results  in  a significant  variation  on  the  magnitude  of 
significant  wave  height  and  parameter  A3  in  shallow 
water.  The  significance  of  this  effect  is  reduced  as  the 
average  water  depth  increases. 

(2)  The  average  zero-crossing  wave  period  slowly 
increases  with  time  during  the  growing  stage  of  the  storm. 
This  is  not  affected  by  the  time  variation  of  water  depth 
regardless  of  the  average  water  depth. 

(3)  The  pattern  of  time  variation  of  parameter 
A 3 (skewness)  during  the  growing  stage  of  a storm  is  always 
very  similar  to  that  of  significant  wave  height  Hs  at 
various  average  water  depths. 

Based  on  these  findings,  it  is  apparent  that  the 
parameter  X3  has  a functional  relationship  with 
significant  wave  height.  The  parameter  \3  computed  from  wave 
records  measured  at  various  water  depths  ranging  from  1.35 
m to  8.77  m is  plotted  against  significant  wave  height  (some 
examples  are  shown  in  Figure  3.4(a)  through  3.4(f)).  The 
results  presented  in  these  figures  are  at  locations  where 
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(a)  Average  water  depth  1.35  m 


(b)  Average  water  depth  3.70  m 


Figure  3.4  Parameter  X 3 as  a function  of  significant 
wave  height 
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(d)  Average  water  depth  5.97  m 


Figure  3.4  (continued) 
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(e)  Average  water  depth  6.97  m 
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Figure  3.4  (continued) 
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waves  were 

measured,  and 

the  water  depth  given 

in 

each 

figure  is 

the  average  water 

depth  during 

the  storm 

at 

that 

location  (see  Figure  A.l 

in  Appendix). 

Included 

in 

each 

figure  is 

the  functional 

relationship 

between 

X 3 

and 

significant 

wave  height  expressed  in  the 

form  of  ^ 3 

= a 

b 

(Hs), 

where  H is  the  significant  wave  height  and  coefficient  "a" 
and  "b"  are  obtained  by  drawing  the  average  line  through 
data  presented  in  the  figure. 

Coefficients  "a"  and  "b"  obtained  from  Figure  3.4(a) 
through  3.4(f)  are  then  plotted  against  the  average  water 
depth.  The  results  are  shown  in  Figures  3.5(a)  and  3.5(b) 
respectively.  These  figures  show  that  the  values  of 
coefficients  "a"  and  "b"  evaluated  at  a location  where  the 
average  water  depth  was  in  the  range  of  5.5  m to  6.0  m 
deviate  from  the  average  line  drawn  in  the  figure.  That  is, 
for  a given  significant  wave  height,  the  skewness  evaluated 
from  the  record  is  greater  than  that  evaluated  by  using  the 
value  from  the  average  line.  This  implies  that  the  water 
depth  of  5.5  m to  6.0  m appears  to  be  the  depth  where  wave 
breaking  occurred  during  the  storm  . This  has  also  been 
mentioned  by  Aranuvachapun  and  Thornton  (1984)  in  their 
study  of  spatial  and  temporal  transformation  of  shallow  wave 
energy.  Hence,  the  non-Gaussian  characteristics  of  the  waves 
in  this  particular  location  are  relatively  pronounced. 

The  coefficients  "a"  and  "b"  are  expressed  as  a 
function  of  the  average  water  depth  h,  as  shown  in  Figure 
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AVERAGE  WATER  DEPTH  IN  METERS 


Figure  3.5  Parameter  "a 
water  depth 


and  "b"  as  a function  of 
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3.5,  and  thereby  the  parameter  A3  can  be  evaluated  in  terms 
of  water  depth  and  sea  severity  as  follows: 

0.59 

0.74  h 

X3  = 1.16  exp( -0 . 42  h)  Hs  (3.2) 

where  h = average  water  depth  in  meters 

Hs  = significant  wave  height  in  meters. 


Next,  the  parameter  A 4 is  evaluated  from  the  records 
and  its  values  are  plotted  against  the  parameter  a3,  as 
shown  in  Figure  3.6.  The  results  in  the  figure  show 
several  interesting  features: 

(1)  There  is  considerable  scatter  in  the  value  of  x4 
for  values  of  A3  less  than  0.2.  These  values  are  obtained 
from  measured  records  taken  where  the  water  depth  is 
relatively  deep,  on  the  order  of  15  to  25  meters.  Although 
the  values  of  A4vary  considerably  for  A3  less  than  0.2, 
the  shape  of  the  non-Gaussian  probability  distribution 
appears  to  be  very  close  to  that  of  the  normal  distribution. 
As  an  example.  Figure  3.7  shows  a comparison  between  the 
non-Gaussian  probability  density  function  with  A3  =0.2  and 
A4  =-  0.15  and  the  normal  probability  density  function  in 
the  standardized  form  . The  comparison  shows  no  significant 
difference  between  these  two  probability  density  functions. 
From  the  results  of  many  similar  comparisons,  it  may  be 
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Figure  3.6  Parameter  A4  as  a function  of 
parameter  A 3 
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Figure  3 . 7 


Comparison  between  normal 
non-Gaussian  distribution 
*4  = -0.15 


probability  distribution  and 
with  \ 3 = 0.20  and 
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safely  concluded  that  the  non-Gaussian  presentation  of 
coastal  wave  profiles  (in  standarized  form)  with  a value  of 
X3  less  than  0.2  can  be  approximated  by  the  Gaussian 

probability  density  function  (in  standardized  form). 

(2)  The  largest  value  of  the  parameter  X3  obtained 

in  the  present  analysis  is  1.26  with  the  parameter  X4  of 
1.98.  This  is  observed  in  a severe  sea  of  significant  wave 
height  2.31  m at  an  average  water  depth  of  1.97  m.  Figure 
3.8  shows  the  histogram  of  the  wave  profile  together  with  a 
portion  of  the  time  series  of  wave  record.  As  can  be  seen, 
the  histogram  of  wave  profile  deviates  from  the  normal 

distribution,  and  the  time  history  of  coastal  waves  show  a 

definite  excess  of  high  crests  and  shallow  troughs. 

(3) There  is  some  scatter  in  the  X4 -values  for  a 


specified 

* 3 * 

For 

example,  the 

value 

of 

X4  spreads  from 

0.45  to 

0.9 

with 

an  average  value  of 

0.7 

for 

X 3 —0 . 8 . 

However, 

the 

difference  in  the 

X4  -value 

for  a 

specified 

X3  does  not  seriously  affect  the  shape  of  the  non-Gaussian 
probability  density  function.  Figure  3.9  shows  a comparison 
between  three  probability  density  fuctions  having  the  same 
X 3 -value ( 0 . 8 ) , but  with  values  of  X4 of  0.45,0.70  and  0.90. 
Included  also  in  the  figure  is  the  Gaussian  probability 
density  function.  The  shape  of  the  non-Gaussian  probability 
density  functions  computed  by  Equation  (3.1)  differs 
significantly  from  that  of  the  Gaussian  probability  density 
function,  but  there  is  no  substantial  difference  between 
the  shapes  of  the  non-Gaussian  distributions.  In  other 
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Figure  3.8  Histogram, :probability  distribution  and  portion  of 
wave  record  for  significant  wave  height  2.31  m at 
average  water  depth  1.97  m 
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Figure  3.9  Comparison  between  normal  probability  distribution  and 
non-Gaussian  distribution  with  X3  = 0.80  and  X4=  0.45, 
0.70,  and  0.90 
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words,  the  scatter  in  the  X4  -value  for  a given  X3 
does  not  result  in  any  serious  difference  in  the  shape  of 
the  probability  density  function. 

From  the  discussion  and  subsequent  conclusion  stated 
in  Item  (1)  and  (3)  above,  the  relationship  between  X4  and 
X3  can  be  formulated  by  taking  the  average  value  of 

X4  for  a specified  X3 -value.  That  is. 


Xu  =< 


1.17 

-0.15  + 1.10  ( X 3 -0 . 2 0 ) for  0.2  $ X3  $0.5 


1.17  1.47 

■0.15+  1.10  ( X 3 -0 . 20 ) + 1.48  ( X 3 -0.50) 


for  X 3 > 0.5 


(3.3) 


In  summary,  waves  at  a given  water  depth  may  be 
considered  as  a Gaussian  random  process  in  mild  sea 
if  the  parameter  X3  evaluated  by  applying  Equation  (3.2)  is 
less  than  0.2.  However,  in  severe  seas  for  which  the 
parameter  X3  is  greater  than  0.2,  waves  must  be  considered 
as  a non-Gaussian  random  process  represented  by  Equation 

(3.1)  consisting  of  terms  containing  the  parameters  X3  , 

2 

X3  , and  X4 . The  parameters  X3  and  X4  can  be 
evaluated  by  Equations  (3.2)and  (3.3)  respectively. 


45 


It  is  of  interest  to  compare  the  non-Gaussian 
properties  of  the  measured  wave  profile  obtained  in  coastal 
waters  with  results  obtained  from  nonlinear  wave  theory.  As 
an  example,  the  skewness  A3  is  computed  from  deviations  from 
the  mean  value  of  nonlinear  waves  generated  by  the  stream 
function  theory  developed  by  Dean  (1974).  A total  of  40 
cases  are  computed  for  different  wave  conditions 
characterized  by  the  relative  water  depth,  h/LQ , and  the 
relative  wave  height,  h/Lq , where  h is  the  water  depth,  H is 
wave  height,  and  Lq  is  wave  length  in  deep  water.  The 
results  of  the  computations  are  shown  in  Figure  3.10  where 
the  dashed  lines  indicating  various  A 3 -values  ranging  from 
0.3  to  1.5  are  generated  by  connecting  values  of  the  same 
skewness.  Included  also  in  the  figure  are  the  measured  data 
for  A3  equal  to  0.50,  0.70,  and  0.90  obtained  from 
analysis  of  coastal  waves.  In  plotting  these  data  obtained 
in  random  seas  in  the  figure,  the  significant  wave  height  is 
considered  as  H,  and  the  wave  length  corresponding  to  the 
average  zero-crossing  wave  period  is  considered  as  Lq . 

As  can  be  seen  in  Figure  3.10,  for  a given  relative 
wave  height,  H/L0,  A3  based  on  the  stream  function  theory 
increases  with  decrease  in  relative  water  depth,  h/L0.  This 
is  true  for  all  cases  where  A3  is  less  than  the  breaking 
index  curve.  However,  the  measured  data  do  not  show  the  same 
trend.  This  may  be  attributed  to  the  difference  in  the 
definition  of  skewness  used  in  the  analysis.  Needless  to 


RELATIVE  WAVE  HEIGHT,  H /L 
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o 


a X,=  0.90 


RELATIVE  WATER  DEPTH,  h/ L0 


Figure  3.10  Comparison  of  skewness  obtained  from  measured 
data  with  that  based  on  stream  function  theory 
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say,  in  the  strict  sense  skewness  should  be  defined  based  on 
a random  process.  Nevertheless,  most  of  the  measured  data 
falls  into  a zone  bounded  by  relative  water  depth  between 
0.06  and  0.20  and  relative  wave  height  between  0.04  and 
0.07. 


3 . 3 Non-Gaussian  Properties  and  Water  Depth 

As  discussed  earlier,  the  non-Gaussian  characteristics 
of  wind-generated  coastal -waves  depend  on  water  depth  as 
well  as  sea  severity.  That  is,  even  though  the  water  depth 
is  shallow,  waves  may  still  be  considered  to  be  a Gaussian 
random  process  if  the  sea  severity  is  mild.  Therefore,  it 
may  be  of  considerable  interest  to  find  the  limiting  sea 
severity  below  which  wind-generated  v/aves  are  considered  to 
be  Gaussian  for  a given  water  depth. 

In  the  discussion  of  Figure  3.7,  it  was  stated  that 
the  non-Gaussian  presentation  of  coastal  wave  profiles  with 
a value  of  *3  less  than  0.2  can  be  approximated  by  the 
Gaussian  probability  distribution,  where  ^3  was  expressed  as 
a function  of  water  depth  and  sea  severity  in  Equation 
(3.2).  Therefore,  it  may  be  well  to  confirm  whether  or  not 
the  significant  wave  height  which  yields  ^3  =0.2  computed 
from  Equation( 3 . 2 ) for  a given  water  depth  represents  the 
limiting  sea  severity  below  which  the  seas  can  be  considered 
to  be  a Gaussian  random  process. 
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For  this,  histograms  of  the  wave  profiles  obtained  at 
various  water  depths  are  examined  to  find  the  limiting  sea 
severity  (significant  wave  height)  below  which  the 
histograms  approximate  a normal  distribution.  As  an  example. 
Figure  3.11  shows  the  histograms  as  well  as  a portion  of  the 
wave  record  obtained  in  a sea  of  significant  wave  height 
1.62  m at  a location  of  water  depth  6.97  m.  This  is  the 
largest  significant  wave  height  below  which  Gaussian 
characteristics  were  observed  during . the  storm  at  this 
location.  The  computed  value  of  the  parameter  x3  (skewness) 
is  0.19.  A similar  analysis  was  done  at  various  depths  and 
the  results  are  shown  in  Figure  3.12.  Included  also  in*  the 
figure  are  the  lines  indicating  various  A.  3 values  computed 
from  Equation  (3.2).  As  can  be  seen,  the  largest 
significant  wave  height  below  which  Gaussian  characteristics 
are  observed  for  water  depths  up  to  15  m agrees  well  with 
the  line  for  x3  =0.2. 

It  is  noted  that  the  largest  significant  wave  height 
observed  at  a water  depth  of  24.4  m was  4.34  m during  the 
storm,  and  that  all  histograms  of  the  wave  profile  obtained 
at  this  location  approximate  a normal  probability 
distribution.  It  thus  appears  that  the  sea  severity  at  this 
location  did  not  reach  a level  during  the  storm  to  produce 
non-Gaussian  characteristics.  Thus,  although  confirmation 
could  not  be  made  for  a water  depth  of  24.4  m,  it  may 
safely  be  concluded  from  Figure  3.12  that  coastal  waves  in 
seas  for  which  the  parameter  A.  3 (skewness)  is  less  than 


tX) 
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Figure  3.11  Histogram,  probability  distribution  and  portion  of 
wave  record  for  significant  wave  height  1.62  m at 
average  water  depth  6.97  m 


SIGNIFICANT  WAVE-  HEIGHT  IN 
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Figure  3.12  Largest  significant  wave  height  below  which 

non-Gaussian  characteristics  are  expressed  as 
a function  of  water  depth 
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0.2  can  be  considered  as  a Gaussian  random  process 
regardless  of  water  depth. 


3.4  Statistical  Characteristics  of  Wave  Period 

Evaluation  of  the  statistical  characteristic  of  wave 
periods  in  finite  water  depth  from  the  information  of  wave 
spectrum  is  extremely  important  for  design  consideration  of 
coastal  structures. 

For  waves  in  deep  waters  which  are  considered  to  be  a 
Gaussian  random  waves,  the  average  zero-crossing  wave  period 
is  evaluated  most  commonly  by  the  following  equation 
developed  by  Rice  (1945).  That  is 


(3.4) 


where 


S (w)  : 


computed  average  zero-crossing  wave  period 

o 

wRS(w)dw  : nth  moment  of  wave  spectrum 


wave  spectrum 


Since  the  waves  in  finite  water  depth  show  non- 
Gaussian  properties,  the  currently  available  method  for 
evaluating  the  period  may  no  longer  be  valid.  Collins 
(1967),  Chakrabarti  and  Snider  (1974),  and  Chakrabarti 
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and  Cooley  (1977)  found  from  the  analysis  of  storm  wave 
records  that  the  average  zero-crossing  wave  periods  computed 
by  Equation  (3.4)  were  less  than  those  measured  from  the 
wave  record.  Goda  (1974)  analyzed  the  wave  data  in  finite 
water  depth  and  evaluated  that  the  average  zero-crossing 
wave  period  computed  from  the  spectrum  by  Equation  (3.4)  was 
20%  less  than  that  measured  from  wave  records. 

To  clarify  the  statistical  characteristics  of  wave 
period  in  finite  water  depths,  an  analysis  is  carried  out  on 
measured  wave  records  obtained  during  a storm  at  water 
depths  ranging  from  1.35  m to  8.77  m.  That  is,  the  average 
zero-crossing  wave  periods  are  evaluated  by  Equation  (3.4) 
and  the  results  are  compared  with  the  average  periods 
obtained  directly  from  wave  records.  The  latter  is  denoted 
bY  Tom 

The  comparisons  of  Toc  and  Tom  made  at  water  depths 
of  1.35  m,  3.70  m,  and  8.77  m in  various  sea  severities  are 
shown  in  Figure  3.13.  Included  also  in  this  figure  are  the 
average  lines  drawn  through  the  data. 

It  can  be  seen  in  the  figure  that, 

(1)  In  general,  the  zero-crossing  wave  periods 
computed  by  Equation  (3.4),  Toc,  are  smaller  than  that 
measured  from  the  wave  record,  Tom.  Toc  appears  to  agree 
with  Tom  only  for  small  periods  (2  sec  or  less)  irrespective 
of  water  depth. 

(2)  The  discrepancy  between  the  magnitude  of  Toc 
and  Tom  is  particularly  pronounced  for  waves  in  shallow 


COMPUTED  AVERAGE  ZERO- CROSSING  PERIOD,  T0  IN  SEC. 
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Figure  3.13  Comparison  between  the  average  zero-crossing 

wave  period  T and  T 

oc  om 
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water.  The  discrepancy  is  reduced  with  increase  in  water 
depth,  and  T appears  to  agree  with  T in  deep  water. 

(3)  Since  the  wave  period  generally  increases  with  the 
increase  of  sea  severity  (significant  wave  height),  it  can 
be  said  that  the  difference  in  magnitude  between  Tqc 
and  T increases  with  the  increase  of  sea  severity  at  a 
given  water  depth. 

Based  on  the  average  lines  shown  in  Figure  3.13,  the 
ratios  of  Toc  and  T^m  are  computed  and  plotted  against  the 
water  depth  as  a function  of  the  computed  value  as 

shown  in  Figure  3.14.  This  figure  may  serve  to  provide 
information  as  to  the  limiting  water  depth  above  which  the 
computed  average  periods  agree  with  the  measured  ones.  For 
example,  the  average  period  of  5 sec.  computed  by  Equation 
(3.4)  can  be  used  for  design  if  the  water  depth  of  the 
location  is  greater  than  9 meters.  Otherwise  the  computed 
value  has  to  be  multiplied  by  the  factor  read  off  from  the 
curve  given  in  the  figure. 

In  summary,  the  average  zero-crossing  wave  period 
evaluated  from  the  information  of  wave  spectrum  by  Equation 
(3.4)  is  smaller  than  that  measured  from  the  wave  record,  in 
general.  The  difference  in  magnitude  between  them  depends  on 
the  water  depth  as  well  as  sea  severity. 
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Figure  3.14  The  ration  Tom/Tocas  a function  of  water  depth 


CHAPTER  IV 


STATISTICAL  PREDICTION  OF  COASTAL  WAVE  HEIGHTS 

In  general,  waves  in  deep  water  are  commonly 

considered  to  be  a Gaussian  random  process  with  a narrow- 
banded  spectrum;  and  hence,  the  statistical  characteristics 
of  wave  amplitudes  and  wave  heights  can  therefore  be 
predicted  through  the  Rayleigh  distribution. 

However,  one  of  the  most  distinct  statistical 
characteristics  of  waves  in  finite  water  depth  is  that  time 
histories  show  a definite  excess  of  high  crests  and  shallow 
troughs  and  cannot  be  considered  a Gaussain  random  process. 
Hence,  the  current  method  for  predicting  wave  heights  based 
the  concept  of  Gaussian  random  process  is  no  longer  valid 
for  probabilistic  prediction  of  wave  heights  in  finite  water 
depth. 

Furthermore,  the  magnitude  of  wave  peaks  in  finite 
water  depth  is  generally  much  greater  than  that  of  troughs. 
Therefore,  probabilistic  characteristics  of  peaks  and 
troughs  for  waves  in  finite  water  depth  should  be  evaluated 
independently  in  contrast  to  the  prediction  method  for  ocean 
waves  in  which  probabilistic  characteristics  of  peaks  and 
troughs  are  assumed  to  be  the  same. 
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This  chapter  presents  a method  to  statistically 
predict  peaks,  troughs,  and  peak-to-trough  excursions  (wave 
heights)  of  waves  in  finite  water  depth.  The  probability 
density  functions  are  developed  based  on  a non— Gaussian 
probability  distribution  with  parameters  which  depend  on 
the  sea  severity  (significant  wave  height)  and  water  depth. 

This  chapter  consists  of  three  sections.  The  first 
two  sections  discuss  the  method  for  predicting  peaks  and 
troughs  as  well  as  peak-to-trough  excursions  of  waves.  The 
third  section  presents  the  application  of  this  prediction 
method  for  estimating  the  significant  value  and  extreme 
value  of  peaks,  troughs  and  peak-to-trough  excursions  of 
waves  in  a given  sea  severity  (significant  wave  height)  at  a 
specified  water  depth.  Numerical  computations  are  carried 
out  for  various  water  depths  ranging  from  1.35  m to  8.77m. 
The  results  are  compared  with  measured  data  obtained  during 
the  ARSLOE  Project. 

4.1  Probability  Distribution  of  Peaks  and  Troughs 

As  stated,  wave  profiles  observed  in  finite  water 
depth  show  a definite  excess  of  high  crests  and  shallow 
troughs.  Hence,  the  probability  density  function  applicable 
to  peaks  and  troughs  should  be  derived  independently.  We 
first  consider  the  probability  density  function  for  the 
peaks . 
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In  this  study,  we  assume  that  the  wave  spectrum  is 
narrow-banded  and  hence  there  exists  only  a single  peak 
during  every  half-cycle  period.  The  probability  that  the 
peak,  which  is  a random  variable  denoted  by  E exceeds  a 
specified  level  E,+  is  given  as  the  expected  value  of  the 
ratio  of  number  of  peaks  exceeding  the  level  £+  to  the 
number  of  peaks  in  the  positive  domain.  That  is 


(4.1) 


where 


F(  £+) 


N 

+ 


cumulative  distribution  function  of  S 

+ 

number  of  peaks  exceeding  a level  E,+  per 
unit  time. 

number  of  peaks  on  the  positive  domain  per 
unit  time, 
expected  value  of  N 
expected  value  of  N 


It  is  noted  that  the  expected  value  E[  N^+/N+  ] is 
equal  to  the  ratio  of  the  individual  expected  values, 
denoted  by  N^+  and  N+  respectively,  only  if  we  assume 
that  the  random  variables  ( N^  + / N+  ) and  N+  are 

statistically  independent.  However,  in  general,  this 

assumption  appears  to  be  acceptable  for  wind-generated 


waves . 
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For  a narrow-banded  random  process,  the  expected 
values  of  and  N+  can  be  approximately  evaluated  from 
the  joint  probability  density  function  of  wave  surface 
elevation  and  velocity  (time  derivative  of  wave  surface 
elevation).  They  are 


N = 
5 + 


0, 


x ) dx 


x ) dx 


(4.2) 


From  (4.1)  and  (4.2),  the  probability  density  function 
of  5 + can  be  derived  as 


d 


d?  + 


Ix|  f ( 5 


+' 


x ) dx 


f(  0, 


x 


) dx 


(4.3) 


In  the  joint  probability  density  function  of  wave 
surface  elevation  and  velocity  shown  in  Equation ( 4 . 3 ) , the 
following  assumptions  are  considered: 

(1)  The  wave  surface  elevation  and  the  velocity  are 
statistically  independent.  To  examine  this  assumption,  a 
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statistical  analysis  on  the  measured  wave  profiles  and 
velocities  obtained  by  differentiating  the  wave  profile  has 
been  carried  out.  In  this  analysis,  the  correlation 
coefficient  of  the  surface  elevations  and  velocities  is 
evaluated . 

Table  4.1  shows  an  example  of  correlation  coefficients 
evaluated  at  three  locations  of  different  water  depths.  The 
table  shows  that  the  magnitude  of  the  correlation 

coefficients  of  x and  x are  very  small  (on  the  order 
of  0.2  or  less)  regardless  of  water  depth  and  sea 

severity.  Hence,  it  may  safely  be  assumed  that  wave  surface 
elevation  x and  velocity  x are  statistically  independent. 
This  results  in  the  joint  probability  density  function  of 
x and  x shown  in  Equation ( 4 . 3 ) being  expressed  as  the 
product  of  the  individual  probability  density  functions  of 
f(x)  and  f(x).  Here,  the  probability  density  function  of 
wave  surface  elevation  f(x)  is  given  in  Equation  (3.1). 
Based  on  the  results  of  analysis  on  the  effect  of  parameters 

on  the  probability  density  function  presented  in  Chapter 

2 

III,  the  terms  with  the  parameters  I3,  A.  3 , and  A4  in  the 
series  are  considered. 

(2)  The  probability  density  function  of  velocity  f(x) 

is  approximately  a Gaussian  distribution  with  zero  mean 

2 

and  variance  . Results  of  a statistical  analysis  of 

velocity  x obtained  from  differentiating  measured  wave 
profiles  have  indicated  that  the  velocity 


approximately 


61 


Table 

4.1  Correlation 
and  velocity 

coefficient 

X 

of  wave 

surface  elevation 

GAGE 

WATER  DEPTH 
(IN  M.) 

SIGNIFICANT 
WAVE  HEIGHT 
(IN  M.) 

SKEWNESS 

pxx 

615 

1.35 

1.58 

0.94 

-0.191 

615 

1.35 

1.83 

1.00 

-0.191 

615 

1.35 

2.05  ' 

1.15 

-0.184 

645 

3.70 

2.10 

0.50 

-0.159 

645 

3.70 

2.55 

0.93 

-0.155 

625 

8.77 

2.93 

0.38 

-0.155 

625 

8.77 

3.75 

0.69 

-0.137 

625 

8.77 

3.63 

0.71 

-0.123 

P xi  : Correlation  coefficient  of  elevation  x and  velocity  x 
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follows  the  Gaussian  distribution.  Examples  of  comparisons 
between  the  histograms  and  the  Gaussian  distribution  of 
velocity  are  shown  in  Figure  4.1.  The  figure  pertains  to 
waves  at  a water  depth  of  3.70  m (12.1  ft.)  in  a sea  of 
significant  wave  height  2.45  m (8.04  ft.). 

For  the  probability  density  function  specifically 
applicable  to  peaks  of  non-Gaussian  random  waves,  it  is 
necessary  to  consider  a variance  which  is  affiliated  only 

with  the  positive  part  (upper  portion)  of  wave  profile. This 

2 

variance,  denoted  by  crx  + can  be  evaluated  by  taking  the 
second  moment  of  the  probability  density  function  of  the 
wave  profile,  f(x),  truncated  at  x=0.  That  is. 


/ 


x f(x)  dx 


a = 

x + 


/ 


f(x)  dx 


(4.4) 


By  applying  the  probability  density  function  given  in 

2 

Equation  (3.1)  to  Equation  (4.4),  the  variance  a x+  is 

expressed  as  a function  of  the  parameter  A3  and  the 

2 

variance  a That  is 

X 
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Figure  4.1  Comparison  of  wave  velocity  histogram  and  Gaussian 
distribution : 

Average  water  depth  3.70  m,  significant  wave  height 
2.45  m 
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Taking  these  conditions  into  consideration,  the 
probability  density  function  of  peaks  can  be  obtained  as 
follows : 


0^  5 < 00 


where 


(4.6) 


*4  c 2 
=1+--  24  *3 


It  is  noted  that  the  probability  density  function 
f ( x ) given  in  Equation  (3.1)  is  not  symmetric  with  respect 
to  x =0,  and  hence  the  differentiation  involved  in  the 
numerator  of  Equation ( 4 . 3 ) results  in  a finite  value  of  f(x) 
0 . To  overcome  this  difficulty,  the  Hermite 


at  S+  = 


polyminal,  H2 


(•H 

' x + 


in  Equation  (4.6),  is  modified  such 


that  it  becomes  zero  for  small  5 


That  is 
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The  probability  density  function  of  troughs  can  be 
derived  through  the  same  procedure  used  for  the  derivation 
of  the  probability  density  function  of  wave  peaks.  The  form 
of  the  probability  density  function  of  troughs  is  like  that 

given  in  Equation  (4.6)  except  5+  should  be  negative  and 

2 

the  variance  <*x  + should  be  replaced  with  the  variance 

2 

°x-  which  is  affiliated  with  the  negative  part  (lower 

2 

portion)  of  wave  profile.  The  variance  ax_  is  given  by 


2 


a 

X - 


f(x)  dx 


(4.8) 


Then,  the  probability  density  function  of  troughs  is  given 


as  follows: 
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where 


-00  < 0 

(4.9) 


L = 


The  Hermite  polyminal  H2  ^ j shown  above,  has  to 
be  modified  in  order  to  let  it  be  zero  for  small  ^ . That 
is 


* 


0 


1 


for  A £ -1 


(4.10) 

for  > -1 

a ' 


Note  that  if  the  parameters  x3  and  X4  are  zero, then 
the  probability  density  functions  of  peaks  and  troughs, 
shown  in  Equations  (4.6)  and  (4.9)  respectively,  become 
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identical  and  reduce  to  the  Rayleigh  probability 
distribution  which  is  applicable  to  peaks  and  troughs  of  a 
Gaussian  random  process. 


4. 2 Probability  Distribution  of  Wave  Heights 

The  probability  distribution  of  the  peak— to— trough 
excursion,  denoted  by  H,  can  be  obtained  as  a distribution 
for  the  sum  of  two  random  variables:  £+  (peak) 
and  £_  (trough).  That  is 


H - 5+  + |CJ 


(4.11) 


where 


H : peak-to-trough  excursion  (wave  height) 
C+  : wave  peak 
£_  : wave  trough 


In  the  derivation  of  the  probability  density  function 
of  wave  heights,  peaks  £+  and  troughs  £_  are  considered 
statistical ly  independent.  This  assumption  appears  to  be 
acceptable  particularly  for  waves  in  shallow  water.  The 
results  of  analysis  of  waves  measured  at  water  depths  of 
1.35  m and  1.97  m have  shown  that  the  correlation 
coefficients  computed  between  peaks  and  troughs  are 
extremely  small,  on  the  order  of  0.15  or  less. 

Thus,  the  joint  probability  density  function  of  wave 
peaks  and  troughs  is  given  as  the  product  of  the 
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probability  density  functions  of  peaks  and  troughs  given 
in  Equations  (4.6)  and  (4.9),  respectively.  The  probability 
density  function  of  peak-to-trough  excursions  can  then  be 
obtained  by  the  following  convolution  integral; 


f(  H ) 


f( 


S+-  H ) dC  + 


(4.12) 


where,  K is  a constant  associated  with  the  normalization  of 
the  probability  density  function, f ( H ) , and  is  given  by 


f CO 

K = / f ( H ) dH 


(4.13) 


By  applying  Equation  (4.12),  the  significant  wave 
height  and  extreme  wave  height  can  now  be  evaluated.  The 
significant  wave  height,  denoted  by  Hs,  is  given  by 

H f(H)  dH  (4.14) 

where  5*  is  the  value  which  satisfies  the  following 
relationship : 

S*  2 

f ( H ) dH  = 


3 


(4.15) 
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The  probable  extreme  wave  height  in  N observed  waves, 
denoted  by  He , which  is  the  modal  value  of  the  n-th  order 
probability  density  function  is  given  as  the  value  which 
satisfies  the  following  relationship. 


H 


ef(H)  dH  = 1 


1 


N 


(4.16) 


It  is  noted  that  the  significant  value  as  well  as 
extreme  value  of  peaks  and  troughs  of  non-Gaussian  waves 
also  can  be  evaluated  simply  by  replacing  the  probability 
density  functions  in  Equation  (4.14)  and  (4.16)  with  that  of 
peaks  and  troughs  shown  in  Equation  (4.6)  and  (4.9). 


4 . 3 Examples  of  Application 

Before  presenting  examples  of  application  of  the 
method  for  predicting  peaks,  troughs  and  peak-to-trough 
excursion  (include  extreme  value)  of  waves  in  finite  water 
depth,  the  computation  procedure  will  be  summarized  as 
follows : 

Suppose  we  want  to  predict  extreme  wave  height  in  a 
sea  of  significant  wave  height,  Hs  (in  meters),  at  a 
location  where  the  water  depth  is  known  as  h (in  meters). 
First,  an  evaluation  is  made  to  identify  if  waves  in  the 
specified  sea  condition  can  be  considered  as  a Gaussian 
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random  process  at  this  location.  This  can  be  done  either  by 
computing  the  skewness  A 3 by  Equation  (3.2)  or  by  referring 
to  the  limiting  line  for  the  Gaussian  random  process  given 
in  Figure  3.11  . If  skewness  is  less  than  or  equal  to  0.2, 
then  the  waves  in  the  sea  can  be  considered  to  be  a Gaussian 
random  process,  and  hence  the  prediction  can  be  made  by 
applying  the  commonly  known  method  which  is  developed  for 
evaluating  wave  height  in  deep  water. 

If  the  computed  value  of  skewness,  ^3  is  greater 
than  0.2,  then  the  waves  have  to  be  considered  as  a non- 
Gaussian  random  process,  and  the  probability  density 
functions  derived  in  the  previous  sections  can  now  be 
applied  for  the  prediction  of  peaks,  troughs  and  wave 
heights . 

In  addition  to  the  two  parameters  ^3  and  ^4  which 

can  be  evaluated  by  Equation  (3.2)  and  (3.3)  respectively, 

2 

the  variance  ax  is  also  necessary  for  computing  the 
probability  density  functions.  Since  the  relationship 

between  the  specified  significant  wave  height  and  the 

2 

variance  o x remains  unknown  in  finite  water  depth,  the 
following  method  is  used; 

2 

(1)  Let  the  relationship  between  H and  variance  0 

2 * 2 

in  finite  water  depth  be  expressed  as  a =(  H /C  ), 

x s 

where  C is  a constant.  As  an  initial  value  of  the 

2 2 

variance,  compute  °x  =(  Hs/4.0  ) which  is  applicable  to 

2 2 

Gaussian  random  process.  Then,  evaluate  a and  a 

X + X- 

by  Equation  (4.5)  and  (4.8),  respectively. 
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(2)  Compute  the  probability  density  functions  of  wave  peaks 

and  troughs,  shown  in  Equation  (4.6)  and  (4.9),  with 

2 2 

the  computed  parameters  ^3  and  ^4  and  crx+  or  ax_ 
respectively. 

(3)  Compute  the  probability  density  function  of  wave 
height  shown  in  Equation  (4. 12), and  then  the  significant 
wave  height  is  computed  by  Equation  (4.14). 

(4)  Compare  the  computed  significant  wave  height  with 

the  specified  value.  If  these  two  values  do  not 
agree,  then  compute  the  variance  again  by  adjusting  the 
constant  C.  Repeat  the  procedure  until  the  computed 
significant  value  agrees  with  the  specified  value. 

(5)  Significant  value  and  extreme  value  of  the  peaks  and 
troughs  are  computed  from  the  final  probability 
density  function. 

As  an  example  of  application,  the  probability  density 
functions  of  peaks,  troughs,  and  peak-to-trough  excursions 
wave  heights  are  computed  at  a water  depth  of  8.77  m in  sea 
severity  of  significant  wave  height  3.39  m.  The  results  of 
computation  are  compared  with  the  measured  data  obtained 
during  ARSLOE  project  as  shown  in  Figure  4.2  through  Figure 
4.4.  Included  also  in  these  figures  is  the  Rayleigh 
probability  distribution  which  is  commonly  used  for  waves  in 
deep  water. 

As  can  be  seen  in  Figure  4.2  and  Figure  4.3,  the 
probability  density  functions  of  peaks  and  troughs  computed 
by  Equations  (4.6)  and  (4.9)  respectively,  agree  reasonably 
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Figure  4.2  Comparison  between  observed  histogram  of  peaks,  Rayleigh 
distribution,  and  newly  developed  probability 
distribution : 

Average  water  depth  8.77  m,  significant  wave  height  3.39  m 


PROBABILITY  DENSITY  IN  I/M. 
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Figure  4.3  Comparison  between  observed  histogram  of  troughs, 
Rayleigh  distribution,  and  the  newly  developed 
probability  distribution: 

Average  water  depth  8.77  m,  significant  wave  height 
3.39  m 
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Figure  4.4  Comparison  between  observed  histogram  of  peak-to-trough 
excursions,  Rayleigh  distribution,  and  the  newly 
developed  probability  distribution: 

Average  water  depth  8.77  m,  significant  wave  height 
■3.39  m 
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well  with  the  histograms  constructed  from  measured  data, 
while  the  Rayleigh  probability  distribution  substantially 
underestimates  the  peak  values  and  overestimates  the  troughs 
value . 

Figure  4.4  shows  that  the  difference  between 
histograms,  Rayleigh  probability  distribution  and  the 
probability  density  function  of  peak-to-trough  excursions  is 
less  pronounced  than  that  observed  for  the  peaks  and 
troughs  independently.  But  the  probability  density  function 
computed  by  Equation  (4.12)  still  shows  a better  agreement 
with  the  histogram  than  that  of  Rayleigh  distribution. 

Based  on  the  probability  density  functions  shown  in 
Figures  4.2,  4.3  and  4.4,  the  significant  value  as  well  as 
extreme  value  can  be  predicted  by  applying  Equation  (4.14) 
and  (4.16).  Comparisons  of  predicted  and  measured 
significant  value  as  well  as  extreme  values  for  peaks, 
troughs,  and  peak-to-trough  excursions  (wave  heights)  is 
tabulated  in  Table  4.2.  It  is  noted  that  the  predictions  are 
made  for  the  same  sea  severity  as  the  measured  one;  hence, 
the  measured  and  predicted  significant  wave  height  are 
identical . 

As  can  be  seen  in  Table  4.2,  the  predicted  significant 
and  extreme  values  for  peaks  and  troughs  computed  by  the 
present  method  agree  reasonably  well  with  the  measured 
values.  However,  there  is  a substantial  difference  between 
measured  and  predicted  values  computed  by  applying  the 
Rayleigh  probability  distribution.  More  comparisons  of  the 
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Table  4.2  Comparison  of  observed  and  predicted  significant 
and  extreme  value  at  average  water  depth  8.77  m 


PEAKS 


TROUGHS 


PEAK-TO-TROUGH 

EXCURSIONS 


SIGNIFICANT  VALUE 


Observed 

2.21 

m 

1.36 

m 

3.39 

m 

Predicted  by 

Present  theory 

2.32 

m 

1.36 

m 

3.39 

m 

Rayleigh 

distribution 

1.70 

m 

1.70 

m 

3.39 

m 

EXTREME  VALUE 


Observed 


3.44  m 


1.73  m 


5.17 


m 


Predicted  by 
Present  theory 


4.23  m 


1.80  m 


5.33  m 


Rayleigh 

distribution 


2.69  m 


2.69  m 


5.37  m 
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predicted  and  measured  significant  values  as  well  as  the 
extreme  values  are  made  at  water  depths  of  1.35  m,  3.70  m 
and  8.77  m in  various  sea  severities,  and  the  results  are 
tabulated  in  Table  4.3  and  4.4  respectively. 

Table  4.3  shows  that  the  significant  values  of  peaks 
and  troughs  predicted  by  applying  the  probability 
distribution  developed  in  the  present  study  agree  very  well 
with  those  measured  from  wave  records.  On  the  other  hand, 
the  significant  values  of  peaks  predicted  by  applying  the 
Rayleigh  distribution  are  substantially  lower  than  those 
measured  from  wave  records,  regardless  of  water  depth  and 
sea  severity.  Inversely,  the  predicted  significant  values  of 
troughs  computed  by  applying  the  Rayleigh  distribution  are 
much  greater  than  those  measured  from  wave  records. 

Table  4.4  shows  that  the  extreme  values  of  peaks  and 
troughs  predicted  by  the  present  method  are  also  in  better 
agreement  with  the  measured  extreme  values  than  those 
predicted  by  applying  the  Rayleigh  distribution.  However, 
the  magnitude  of  extreme  peak-to-trough  excursions  predicted 
by  applying  the  present  theory  is  much  greater  than  that  of 
the  observed  excursions  for  waves  in  shallow  water  areas, 
say  less  than  4 meters.  This  may  be  attributed  to  wave 
breaking  phenomenon  frequently  observed  in  shallow  water 
areas.  Since  the  present  theory  does  not  consider  the 
breaking  phenomenon,  the  computed  extreme  values  of  peak-to- 
trough  excursions  may  be  high  for  waves  in  shallow  water 
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Comparison  results  make  it  clear  that  the  application 
of  Rayleigh  probability  distribution  for  peaks  of  non- 
Gaussian  random  waves  in  finite  water  depth  results  in  a 
substantial  underestimation  in  the  prediction  of  significant 
value  as  well  as  extreme  value.  This  underestimation  may  be 
a serious  problem  for  the  design  and  safety  operation  of 
coastal  structures,  especially  for  those  which  are  designed 
to  preclude  wave  overtopping.  In  short,  it  is  highly 
desirable  that  the  probabilistic  prediction  of  peaks  and 
troughs  of  waves  in  finite  water  depth  be  based  on  the 
concept  of  the  non— Gaussian  random  process. 


CHAPTER  V 


CONCLUSIONS 

This  paper  presents  the  results  of  a study  on  non- 
Gaussian  characteristics  of  waves  in  finite  water  depths. The 
objectives  of  the  study  are  (1)  to  clarify  the  non-Gaussian 
properties  of  coastal  waves  which  depend  on  water  depth  as 


well  as 

sea  severity. 

and  (2)  to  develop 

a method 

to 

predict 

peaks,  troughs 

and  peak-to-trough 

excursions 

of 

waves . 

To 

achieve  the 

first  objective. 

a statistical 

analysis 

is  carried  out 

on  wave  time  histories  obtained 

at 

various 

water  depths 

during  a storm  by 

the  Coastal 

Engineering  Research  Center  (CERC)  in  the  Atlantic  Ocean 
Remote  Sensing  Land-Ocean  Experiment  (ARSLOE)  Project.  From 
the  results  of  statistical  analysis,  the  following 
conclusions  are  drawn: 

(1)  Time  histories  of  coastal  waves  show  a definite 
excess  of  high  crests  and  shallow  troughs,  in  general,  and 
waves  are  considered  to  be  a non-Gaussian  random  process. 
The  non-Gaussian  properties  depend  on  water  depth  and  sea 
severity. 

(2)  The  results  of  analysis  of  data  show  that  the  non- 
Gaussian  characteristics  of  wave  profiles  can  adequately 
be  represented  by  the  probability  density  function  developed 
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by  Longuet-Higgins  which  is  expressed  in  a series  form  with 
parameters  ^3»  ^4»  ^5/  etc.  It  is  found  in  this  study 

that  the  parameter  which  represents  the  skewness  of  the 
wave  profile  is  the  dominant  parameter  affecting  the  non- 

Gaussian  characteristics  of  coastal  waves,  and  that  the 

2 

combination  of  terms  with  the  parameters  ^3  ^3  and 

A4  (where  ^4  = kurtosis  -3)  in  the  probability  density 
function  best  represents  coastal  wave  profiles. 

(3 ) In  order  to  obtain  the  non-Gaussian  properties 
as  a function  of  water  depth  and  sea  severity,  the  parameter 

^ 3 is  expressed  as  a function  of  water  depth  and  sea 
severity.  Furthermore,  the  functional  relationship  between 
the  parameters  ^3  and  *4  is  derived. 

(4)  The  coastal  waves  may  still  be  considered  to  be  a 
Gaussian  random  process  if  the  sea  severity  is  mild.  It  is 
found  that  the  shape  of  the  non-Gaussian  probability  density 
function  with  the  value  of  *3  less  than  0.2  is  very  close 
to  that  of  a normal  distribution.  Hence,  the  coastal  waves 
for  which  the  parameter  ^3  (skewness)  is  less  than  0.2  can 
be  considered  as  a Gaussian  random  process  regardless  of 
water  depth.  Based  on  this,  it  is  possible  to  evaluate  the 
limiting  sea  severity  below  which  the  coastal  waves  are 
considered  to  be  Gaussian  for  a given  water  depth. 

(5)  The  average  zero-crossing  periods  of  coastal  waves 
evaluated  from  the  most  commonly  used  formula  ( developed 
based  on  a Gaussian  random  process  assumption)  are  smaller 
than  those  measured  directly  from  wave  records. 


This 
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difference  in  magnitude  between  them  is  particularly 
pronounced  in  shallow  water  and  severe  seas. 

For  the  second  objective  of  this  study,  the 
probability  density  functions  applicable  to  peaks,  troughs 
and  peak-to-trough  excursions  of  waves  in  finite  water  depth 
are  developed  based  on  the  findings  obtained  in  the  first 
phase  of  the  study.  The  results  of  numerical  computations 
carried  out  at  various  water  depths  for  various  sea 
severities  show  that  the  probability  density  functions  of 
peaks  and  troughs  both  reasonably  agree  with  the  histograms 
constructed  from  measured  data.  The  Rayleigh  probability 
density  function  which  is  commonly  used  for  predicting 
amplitudes  and  heights  of  waves  in  deep  water  substantially 
underestimates  the  magnitude  of  peaks  and  overestimates  the 
magnitude  of  troughs  of  coastal  waves. 


APPENDIX 


WAVE  DATA 

The  wave  data  used  in  the  present  study  were  obtained 
by  the  Coastal  Engineering  Research  Center  (CERC),  U.S. 
Army,  at  its  Field  Research  Facility,  located  at  Duck,  North 
Carolina.  The  Field  Research  Facility  has  a 550-meter-long 
research  pier  extending  into  the  Atlantic  Ocean,  equipped 
with  seven  Baylor  resistance-type  wave  gages  (see  Figure 
A.l).  Extensive  wave  measurements  were  made  at  this  site  for 
two  months  in  1980  under  the  Atlantic  Ocean  Remote  Sensing 
Land-Ocean  Experiment  (ARSLOE ) Project. 

During  October  23-25,  an  extratropical  cyclon  moved 
directly  through  the  experimental  area  with  wind  speeds  on 
the  order  of  10-15  m/sec.  Continuous  wave  records  were  taken 
by  the  CERC  Field  Research  Facility  during  the  storm.  The 
wave  data  analyzed  in  the  present  study  are  those  recorded 
at  the  seven  locations  along  the  research  pier  and  an 
additional  two  buoys  (wave  riders)  located  offshore  along 
the  extended  line  of  the  pier.  The  location  of  the  wave 
gages  and  wave  riders  as  well  as  the  average  water  depth 
measured  on  October  21  and  27  by  CERC  are  given  in  Figure 
A.l. 
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Figure  A.l  Coastal  Engineering  Research  Center  (CERC) 

Field  Research  Facility,  Duck,  North  Carolina 
(By  courtesy  of  CERC) 
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